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1. Introduction 

Five years ago, an extensive program for the study of strongly interacting matter 
through high energy heavy ion coUisions was initiated at CERN and BNL [1-5] A crucial 
element in this endeavour is the production of systems which arc indeed thcrnializcd - 
and not just an uncorrelatcd superposition of nucleon-nucleon collisions. Thermalisation 
can only be achieved through sufficient rescattering of the incoming nucleons among each 
other and with secondaries, and of the secondaries among each other. It is not obvious 
if, how and at what stage thermalisation is achieved. If it is achieved at all, and if 
it is not destroyed subsequently by possible non-equilibrium features [6-11], then the 
produced hadrons should be observed in ratios determined by the thermodynamics of 
a system at freeze-out temperature T/ and freeze-out baryon density ub- The aim of 
this work is to determine and tabulate these ratios for the various hadron species as 
function of T/ and ub- Our tables should allow a thermalisation analysis of any present 
and future heavy ion data, and we shall illustrate how such an analysis is carried out. 
We should emphasize that we are not providing a model of any kind for the space-time 
development of the systems produced, but merely a check of the degree of thermalisation 
indicated in the observed hadron ratios. In particular, the collisions between light ion 
projectiles and heavy ion targets studied studied up to now very likely do not lead 
to complete thermalisation, and so we expect experiments to deviate from the ratios 
listed here. These deviations should, however, lead to values between the ratios from 
p-p collisions and the thermal ones. If they don't, this could be a hint for some new 
effect. Furthermore, if there is thermalisation, one cannot say anything about how 
it was established. Thermal systems have no memory, and hence the observation of 
thermal hadron production ratios neither excludes nor confirms quark deconfinement at 
an earlier stage. 

It should be noted that the thermalisation we are considering here refers to the 
approach to a statistical equilibrium among the different hadron species. It is conceiv- 
able [11] that such a "chemical" equilibrium sets in later than a "thermal" equilibrium 
among different hadrons of the same species. So far, however, there seem to be no 
clear-cut tests for the latter. 

2. Thermal Hadron Production 

Hadronic matter in equilibrium is specified by three parameters : the temperature 
Tf, the baryon number density ub, and the strangeness density ns; the latter is zero 
unless one is interested in globs of strange matter [10]. Hadronic matter in thermal and 
chemical equilibrium will have a certain composition of particle species, characterized 
by "chemical" potentials. In addition to the baryon number chemical potential fiB, 
we have a strangeness chemical potential //g, since particles of a given strangeness can 
have different baryon numbers (e.g., and A). The value of is fixed by giving 
the overall baryon number density ub, and that of by fixing the overall strangeness. 
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which in our case will be zero. The number densities for all hadron states are then given 
in terms of only hb and Tf. There are many particle species, and if they are produced 
in equilibrium, all their ratios must be fixed in terms of these two parameters. It is thus 
straightforward to test if there was equilibrium, and since this question is so basic, it 
should be one of the first tests to be made. 

We include only the essential interaction features of a multi-hadron system near 
freeze-out: resonance production in general and particle repulsion at high baryon den- 
sity. The bulk of meson production is known to take place through intermediate reso- 
nance states, and we therefore start with an ideal gas containing the observed (strange 
and non-strange) mesonic and baryonic resonances. A crucial question here is how 
many of the observed resonances we should include. There is no clear-cut answer - 
by including more and more resonance states, wc evidently go back further and fur- 
ther in the evolution of the system. To obtain a definite scheme, we shall include all 
resonances up to a mass of 2 GeV. In order to test the importance of high mass contri- 
butions, we have also carried out calculations including only the resonances of the basic 
multiplets, i.e., the pseudoscalar meson octet plus singlet, the vector meson octet plus 
singlet, the baryon octet and the baryon decuplet. - To obtain the distributions of the 
actually observed hadrons, in particular pions and kaons, we then let the resonances 
decay according to the measured branching ratios. This procedure is well known; see 
e.g. [12]. 

For small or vanishing baryon number, this would suffice to fix our model; at high 
ub, however, a strong, short-range repulsion sets in between baryons, even if they 
have diff'erent quantum numbers; this can be described in terms of excluded volumes, 
repulsive interaction contributions, or Pauli blocking. We shall here adopt a hard- 
sphere picture [13], in which the baryons in the system eff'ectively remove a portion of 
the available spatial volume. As a consequence, the baryon number density ub becomes 

= [1 + V,n%] ' W 

where denotes the density calculated for an ideal gas of pointlike baryons. We note 
that for — > oo, ^ 1/^, i.e., we have dense packing* of hard-sphere baryons 
with an intrinsic volume V. In the case of several baryon species a = 1,2, ...,r, eq.(l) 
becomes 

where the baryon volume could be proportional to the mass, as suggested e.g. by 
the bag model. Here we take one universal volume V = 4nR'^/3, with R = 0.8 fm as 
an indicative radius for the nucleon, and apply the excluded volume correction (1) to 
all particle species, mesons as well as baryons. Note that this does not imply repulsion 



* Actually 1/V is the limiting value for incompressible but deformable baryons; a 
true hard-sphere system yields 1/cV, with c > 1 [14]. 
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between mesons; it only reduces space accessible to them, because some is already 
"used up" by the hard-core baryons. For matter at ub = 0, the full volume is once 
more available; for T = 0, we recover the cold nucleon gas with hard-sphere repulsion 
[13]. The correction term moreover cancels out in all particle ratios. - This completes 
our picture. 

To illustrate the resulting pattern, we note that 

_ / _fP I (n^ 

""^^ ~ J (27r)3 [eiEKT^s)/T _ 1] 

gives the number density of directly produced thermal kaons, and 

_ f d^P I / .^ 

~ J (2^ [e(£:A-(MB-/.s))/T + i] ^^j 

that of directly produced thermal, pointlike A's; here Ek — (p^ + m|^)^/^ and Ea = 
(p^ -I- m^)^/^. The actual production rates are obtained from this by including baryon 
repulsion according to eq.(2) and by adding to eqs.(3) and (4) the corresponding val- 
ues from resonance production, such as K^* — > K'^ + 7r~ , using the experimentally 
determined branching ratios. 

Requiring the difference between all strange particle number densities to vanish, 

E^s - = 0' (5) 

i i 

fixes Us and thus gives us all particle densities in terms of T and b ■ Summing at fixed 
fj,B over all baryon contributions, 

E^'s-E^'s = , (6) 

i i 

determines the overall baryon number density ub- Note that in eq.(6) we sum over the 
physicalhaxyon densities as defined by eq.(2). Eqs.(5) and (6) provide a relation between 
the strange chemical potential and the baryon chemical potential; it is shown in Fig.l. 
With increasing temperature, ns is seen to increase faster with fiB, and for temperatures 
T> 200 MeV, Hs/iJ'B is found to be rather constant in the range considered here. In 
particular we note that for T ~ 200 MeV and up to about — 500 MeV, we have 

~ (l/3)/U5. In terms of the baryon number and strangeness chemical potentials for 
quarks, /x^ and /x^, respectively, the hadronic chemical potentials become Hb = '^IJ^b and 
l^s = fJ-b ~ fJ-s- Hence ^jlb — ^l^-s just translates into ~ 0. Since this arises here quite 
naturally in a hadronic picture, it is difficult to consider it as evidence for a primordial 
quark-gluon plasma [15]. 

In Fig. 2 we illustrate the effect of hard-core repulsion between baryons, by plotting 

as function of hb for a hard-core volume V — and for one with radius R — 0.8 
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fm. We note that the repulsion begins to become crucial around ub— 0.1 fm~^, i.e., 
already below standard nuclear density . 

We should emphasize that the particle ratios which we will discuss in this paper 
are by construction independent of the specific relation between the baryon chemical 
potential and the baryon density: as already mentioned, the repulsion introduced in 
equation (1) cancels out in the ratios. The relation between baryon number density and 
baryon chemical potential, on the other hand, does depend on the specific form (1). In 
the tables given at the end of this paper, the baryon number density is determined by 
this form with the hard core volume as defined there. 

3. Thermal Hadron Rates 

In this section we will discuss the general features of the various hadron species. 
Since we are here only interested in the qualitative behaviour, the results for this section 
were obtained using only the basic hadron multiplets, unless otherwise noted. 

3.1 Pions 

Pions are present in thermal equilibrium in the heat bath and in addition in the 
decay products of almost all heavier particles. As long as the temperature is below 100 
MeV, only very few heavy particles are present and almost all pions are of direct thermal 
origin. It is well known, however, that the number of directly produced thermal pions 

decreases at higher temperatures (see e.g. [6, 7], [12], [16-19]). This becomes evident in 
Fig. 3, where the ratio of all pions (i.e., including decay products) to direct thermal pions 
is plotted as a function of temperature. It increases rapidly above 100 MeV. For higher 
temperatures, the total positive pion density is thus given by 

n7r+ =n^+ (thermal) 

+0.285 

2 

+ -np + nu) + 0.134 
1 1 

1 1 1 

+ 2?^A + g^s* + -TIE.* 

1 1 1 

+ 3^A+ g'^E* + g'^H* 

+ ■•-, (7) 

in which the different factors in front of the resonance densities are obtained from the 
corresponding branching ratios into positive pions. As the temperature is increased 
still further, more and more terms will start contributing in the sum in equation (7). 
At fixed temperature, the ratio of all pions to thermal pions furthermore increases 
rapidly as a function of baryon chemical potential. This is shown in Fig. 4 for two 
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different values of the temperature T. The increase is due to the fact that the number 
of baryonic resonances increases with the baryon chemical potential, and these provide 
more pions as decay products. Wc notice that the fraction of thermal pions quickly 
becomes very small; at high temperatures and densities, it is negligible compared to the 
number of pions from decay products of heavier resonances. It is therefore important 
to specify how many such resonances are taken into account in the hadron gas, since 
the number of pions keeps increasing as more and more resonances are included. This 
dependence has been investigated and was found to be important [12,16,20]. In principle 
one must include all resonances; in practice, the contributions from heavy states become 
important only when the temperature of the system is relatively high. We show in Fig. 3 
both the result obtained from our tables, which include all observed resonances up to a 
mass of 2 GeV, and the same ratio obtained by including only the fundamental hadron 
multiplets. We see that up to temperatures of about 150 MeV, including more heavy 
resonances has essentially no effect; at 200 MeV, however, the ratio of all to thermal 
pions has already increased by about 50%. In any case, most pions in the final state of 
a relativistic heavy ion collision are certainly decay products of resonances. The pion 
spectra therefore do tell us something about equilibrium, but only in a very indirect 
way. 

3.2 Kaons 

Kaons are at first sight similar to pions. We see in Fig. 5 that the ratio between 
the total number of kaons and those that are solely of thermal origin also rises rapidly 
above a temperature of 100 MeV, just as is the case with pions, even though the rise 
is not nearly as steep. Decay kaons come mainly from strange mesonic resonances; 
baryons only rarely have kaons in their decay products. Even strange baryons do not 
decay copiously into kaons, witness the baryons in the decuplet, where only the 0~ 
gives kaons. As function of the baryon number density, the ratio of all kaons to thermal 
kaons therefore does not behave like the corresponding ratio for pions; it decreases with 
increasing //b, as seen in Fig. 6. 

The total density of positive kaons is determined by the relation 

nK+ = nK+ (thermal) + 0.495 + 0.5 uk* + • • ■ • (8) 

When the temperature is increased much beyond 150 MeV, more terms would have 
to be added to the right-hand side of the above equation. As mentioned, the baryons 
give almost no contributions, in contrast to the pion case (7), where they contribute 
very much. Nevertheless, we conclude again that a considerable part of the kaons in 
the final state of a relativistic heavy ion collision are not of direct thermal origin, but 
decay products of strange mesonic resonances. The situation here is somewhat less 
pronounced than it was for pions, however, where those of direct thermal origin could 
almost be neglected. The tables at the end of this paper were calculated using equation 
(8) without any further resonance contributions. 
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3.3 Vector Mesons 

Vector mesons have attracted much attention, with the observation of J /if: sup- 
pression compared to the Drell-Yan continuum and of the enhancement of production 
compared to that of p and u [5]. J/t/^ production hes outside present considerations, 
since the J/ifj is too heavy to be produced by thermal colhsions. We therefore turn our 
attention to the ratio 0/ (p + oj), which is shown in Fig. 7 as a function of temperature. 
This ratio is independent of the baryon chemical potential, since non-strange mesons are 
not sensitive to the baryonic content of the hadronic gas. Strange mesons are sensitive 
to riB through the strange chemical potential, e.g. via associated production of A and 
. The (j)-m.esoii is of particular interest, because its quark structure never enters the 
hadronic gas calculations; from a hadronic point of view, the is just a heavy u. Only 
on the quark level does the ss-content enter. - With increasing hb-, the (j)/{p + uj) ratio 
actually drops, because the relative weight of the 4> decreases, as shown in Fig. 8. 

The order of magnitude of the ratio (f)/{p + u>) is given by 

0/(p + w) - ^exp[(mp -m^)/T]. (9) 

For T = 200 MeV, this gives a value of 0.14; the actual distributions lead to slightly 
larger values than the estimate (9). We have not taken into account vector mesons 
originating from the decay of heavy resonances, since the relevant branching ratios are 
not very well established. Their inclusion would lower the ratio, since it will mainly 
enhance the number of p- mesons. 

3.4 Strange Baryons 

Since it is experimentally difficult to separate E° and A, we shall consider these 
together and denote by nx the overall production rate of the two species. Again there 
are non-negligible decay contributions, from the decuplet baryons S* (1.385) and Q~ 
(1.672), which decay into A's and S's. We therefore get 

1 

nx — riA (thermal) -I — ns (thermal) -|- 0.9 n-^* + 0.68 uq + ■ ■ ■ , (10) 

3 

where counts all charge states, so that nj:/3 gives us the neutrals. - In a similar 
way, S gets contributions from its decuplet partner the S* (1.530), leading to 

ns = nE:(thermal) -|- n-,* + • ■ • . (11) 

In Fig.9 we show the relative importance of the S (1.32 GeV) in comparison to the S* 
(1.53 GeV), as a function of temperature for a fixed fiB = 200 MeV. Here we plot the 
ratio (S-|-S*)/S, i.e., (all S/thermal S). Again one observes a quick increase of this ratio 
above temperatures of 100 MeV, similar to the one observed in Fig. 5 for kaons. This 
shows that even a heavy particle like the S is strongly influenced by heavier resonances. 
There thus seems to be no case where heavy resonances can be neglected. Contrary to 
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the kaon case, the ratio considered here is independent of the baryon chemical potential, 
because both particles are baryons with identical quantum numbers and all dependence 
on chemical potentials therefore cancels out. The tables at the end of this paper were 
calculated using all resonances having masses up to about 2 GeV. 

To obtain an order of magnitude estimate of the A/ A ratio, one can use 

(A/A) ~ exp[(-2//B + 2//5)/T], (12) 

while for the S/S ratio one has 

(S/S) - exp[(-2/XB + 4/xs)/T]. (13) 

These estimates neglect the dependence on transverse momentum; in this case, the A 
and the A* (and similarly the S and the S*, as well as the corresponding antibaryons) 
have the same dependence on /ib- Hence all "all/thermal" ratios are here approximately 
^^-independent. However, since the momenta were neglected, these results should be 
used only as rough approximations. Nevertheless, the inequality 

(S/S) > (A/A) (14) 

will hold even for the complete calculation, since the ratio (12) contains a factor 2//^ 
less in the exponent than the ratio (13). 

4. The Analysis of Equilibrium Behaviour 

We will now illustrate for some specific cases how our tables can be used to study 
the degree of equilibration found in heavy ion data, and what conclusions one might 
draw from the results. 

4.1 Strange Baryons 

Here we shall consider the data on strange baryons provided by the WA85 [2] and 
NA35 [3] collaborations at CERN. From two measured ratios we can estimate the two 
thermal parameters, i.e. the temperature and the baryon chemical potential. The data 
of [2] provide the ratios A/A = 0.13 ± 0.03 and S/S = 0.39 ± 0.07. Each leads to a 
"band" in the Tf — hb plane, which can be determined from our tables. The result 
is shown in Fig. 10, where we see that the two bands cross in the rather small shaded 
region. From this we would get the ranges 170 MeV < Tf < 220 MeV for the freeze- 
out temperature and 220 MeV < fiB ^ 420 MeV for the baryon chemical potential, if 
there was equilibrium at the freeze-out of these baryons. Note that the temperature 
range is roughly in accord with temperature estimates from the transverse momentum 
distributions, which from large transverse momentum data give Tpt ~ 230 — 240 MeV 
[2]. We should further emphasize that the pT cut of the experiment does not matter 
here, since in both ratios the masses are the same in numerator and denominator. 
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So far we have used two data points to estimate two parameters. To see if these 
really correspond to the freeze-out parameters of an equilibrium hadron gas, we need 
further data. These are provided by results for S~/A and E]~/A [2]. Due to the 
different masses involved, these ratios have to be corrected for the px cuts made in the 
experiment. This reduces the measured ratios, since the S-particles have a heavier mass 
and thus start at higher values of the transverse mass tut = (py -|-m^)^/^. We estimate 
the corresponding correction factor to be 

^ rriTdmT exp{-mT/T) ^^^^ 
(mT-cut) XTa ^TdrriT exp{-mT/T) 

For equal masses, this correction factor is unity; this is the case for the S/S and the 
A/A ratios. For S~/A and S~/A, the correction factor is approximately 0.48. Using 
the slopes given by the WA85 collaboration, this leads to S~/A = 0.10 ± 0.02 and 
S /A = 0.29 ± 0.10 for the corrected ratios. 

As above, we now calculate the band in the T — hb plane corresponding to E]~/A = 
0.10 ± 0.02; in Fig. 11 it is combined with the results shown in Fig. 10. We see that 
this ratio is also compatible with the T — hb range obtained in Fig. 10. For S~/A = 
0.29 ±0.10, we include in in Fig. 11 only the line corresponding to the value 0.29, without 
errors, in order to keep the graphics clear. It is seen to be compatible with the region 
of Fig. 10 as well. We therefore conclude that the strange baryon ratios as measured 
by [2] are quite consistent with equilibrium at freeze-out for T = 200 ± 20 MeV and 
IJ,B = 320 ± 100 MeV. Note that this value of fiB corresponds to ub— 0.3 fm~^, i.e., to 
about twice nuclear density. 

Let us now turn to the ratio A./K° = 0.76 ± 0.16, measured by [3]. The resulting 
line is included in Fig. 11; it is also absolutely consistent with the previously obtained 
freeze-out parameters, even though we are now looking at the ratio of a strange baryon 
to a strange meson. 

The picture becomes less consistent, however, if we compare the A and A to all 
negative hadrons, i.e., to mainly non-strange mesons. The reported measurements [3] 
A/h~ = 0.08 ± 0.01 and A/h~ = 0.015 ± 0.005 lead to the cross-over region shown in 
Fig. 12. Although the overlap regions in Figs. 11 and 12 just touch, the region in Fig. 12 
does point to lower values of the freeze-out temperature as well as baryonic chemical 
potential. However, before drawing any conclusions, let us consider what the ratios 
among only mesons indicate. 

4.2 Strange and Non-Strange Mesons 

The ratio (f)/{p+uj) was determined in [5] and [21]; both get a value around 0.3 - 0.4 
at the highest energy density, and as seen in the dimuon channel. Correcting this by the 
dimuon branching ratios for (0.025), p and uj (0.007), we get from [21] the value 0.11 
± 0.01 for the ratio 4>/{p + Lj). Since all three mesons are non-strange, the thermal ratio 
is independent of hb- In Fig. 7 we had shown the behaviour of the thermal 0/(p + u>) 
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as function of Tf. The measured value is included in the figure; is seen to require very- 
low temperatures, between 100 and 120 MeV. It is certainly not compatible with the 
Tj~200 MeV found above in the strange baryon analysis. 

This conclusion is supported by an analysis of K/n ratio. In [3], the ratio K"/n~, 
equivalent to {K^ + K~)/ (tt"*" + tt"), is found to be 0.11 ± 0.02 for the integrated data, 
with a peak value of 0.15 ±0.04 at midrapidity. The resulting bands in the Tf—^B plane 
are shown in Fig. 13; with values between 100 and 120 MeV (for the integrated data), 
they also require a considerably lower temperature than found from strange baryons. 
The same conclusion had been reached in an analysis of the BNL data [4] for the K'^ /tt'^ 
ratio [22-24], where temperatures around 100 MeV were needed. 

4.3 Conclusions 

We have thus found two sets of data, each in itself compatible with equilibrium 
hadron ratios, but leading to different freeze-out parameters. Before we comment on 
possible interpretations, let us note some caveats. We have compared experiments cover- 
ing different kinematic regions, and we have considered results integrated over kinematic 
regions showing quite non-uniform behaviour (mid-rapidity vs. forward/backward, for 
example). This is not really meaningful, and hence everything done here should only 
serve as an illustration. Ideally, one should use ratios from symmetric nuclear collisions 
(A — A) at mid-rapidity and integrated over transverse momenta from zero to 1.5 - 2 
GeV. Hopefully, the ^tt-beam at BNL and the P6-beam at CERN will provide such 
data. - The measurements of vector mesons through dilepton or photon decays may 
lead to incorrect rate estimates, if the mesons are short-lived enough to decay in the 
medium [25-28]. 

If we take the result of our analysis serious in spite of these caveats, they would 
seem to indicate that the observed strange baryons freeze-out at an earlier and hence 
"hotter and denser" stage than the mesons. It would be interesting to check if this is 
also reflected in different temperatures from transverse momentum spectra; there seem 
to be hints in that direction. 

In summary: in this paper we have determined and tabulated the ratios for the 
various hadron species as function of the freeze-out temperature Tf and the baryon 
number density n^, or equivalently, the baryonic chemical potential fiB- These tables 
should allow a study of the degree and temporal succession of thermalisation for the 
observed hadron species. 
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Figure Captions 



Fig. 1: Strange chemical potential //g as a function of baryon chemical potential hb for 
temperature values from T = 150 MeV to T = 300 MeV. The solid line was 
calculated for a temperature of 10 GeV. 

Fig. 2.: Baryon density ub as function of baryon chemical potential /Us at temperature 
T = 200 MeV for pointlike baryons (dash-dotted line) and for baryons with hard- 
core repulsion (full line). 

Fig. 3: Ratio of all pions to thermal pions as function of temperature T for baryon chem- 
ical potential iib = 200 MeV. The full line includes the contributions from the 
fundamental hadron multiplets only, the dash-dotted line contributions from all 
resonances with masses up to 2 GeV. 

Fig. 4: Ratio of all pions to thermal pions as function of baryon chemical potential jiB 
for temperatures T = 170 MeV (fuU line) and T = 200 MeV (dash-dotted line); 
calculations include the fundamental hadron multiplets only. 

Fig. 5: Ratio of all kaons to thermal kaons as function of temperature T for baryon chemical 
potential fJ,B = 200 MeV; calculations include the fundamental hadron multiplets 
only. 

Fig. 6: Ratio of all kaons to thermal kaons as function of baryon chemical potential fiB 
for temperatures T = 170 MeV (full line) and T = 200 MeV (dash-dotted line). 
Calculations include the fundamental hadron multiplets only. 

Fig. 7: Ratio of cf) to p+u) as function of temperature T; calculations include all resonances 
up to masses of 2 GeV. 

Fig. 8: Ratio of (/> to p -|- a; as function of baryon chemical potential /x^ at temperature 
T = 200 MeV, including contributions from all resonances with masses up to 2 
GeV. 

Fig. 9: Ratio of all S's to thermal S's as function of temperature T for baryon chemical 
potential (Ib = 200 MeV; calculations include the fundamental hadron multiplets 
only. 

Fig. 10: T — jiB region consistent with the measured ratios for A/A and S/S; the experi- 
mental error limits of the ratios are indicated. 

Fig. 11: T — fiB region consistent with the measured ratios for A/A, S/S, E]~/A, S /A, 
and A/K°; the experimental error limits of the ratios are indicated, except for the 
last two ratios. 

Fig. 12: T — fiB region consistent with the measured ratios for A/h~^ (dashed line) and for 
A/h~ ratio (full line); the experimental error limits of the ratios are indicated. 

Fig. 13: T — fiB region consistent with the measured ratio for Kg/h~; the experimental 
error limits of the ratio are indicated. 
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